Linking Classical and Quantum Stochastic Processes 



Vahid Karimipour and Lalch Mcmarzadeh 
Physics department, Sharif University of Technology, P.O. Box 11155-9161, Tehran, Iran 

We define a map which relates four dimensional classical stochastic matrices to qubit quantum 
channels. The map preserves the spectrum and the composition of processes. To do this we introduce 
the concept of Bloch tetrahedron which plays the same role in the classical context as the Bloch 
sphere in the quantum context. A similar map is also induced between dynamical generators of 
classical and quantum stochastic processes. Possibilities for generalization to arbitrary dimensions 
are also discussed. 
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The state of an open physical system interacting with 
an environment on which we have little or no control, is 
given by a probability vector \P) in the classical domain 
or a density matrix p in quantum domain. Linear evo- 
lution equations describe the dynamics of such states in 
successive time steps. In the classical case, the dynamics 
is given by \Pn+i) = Q\Pn) [H where n denotes the time 
step and the stochastic matrix Q collects the jumps prob- 
abilities. In the quantum domain, the dynamics is repre- 
sented by a completely positive trace-preserving (CPT) 
map or quantum channel Pn+i = £iPn) HI- These equa- 
tions provide an abstract but very general framework in 
which the dynamics of many different phenomena both 
in classical and quantum physics can be analyzed. 

In classical domain, this might describe a random walk 
by a particle, diffusion and reaction of particles, succes- 
sive conformation of a polymer chain, radioactive decay 
of nuclei, or population dynamics of any system of inter- 
acting entities. Its application even goes beyond physics 
to other fields of science, like biology and finance [3[. 

In quantum domain such an equation may describe an 
application of a gate in a quantum processor, a communi- 
cation channel in time (quantum storage device) or space 
(an optical fiber), or any open systems dynamics where 
one merely has partial access to the relevant degrees of 
freedom 

How quantum stochastic processes are related to 
classical ones? Is it possible to construct a quantum 
channel for every stochastic matrix? How about the 
converse? If yes, how the properties of the classical 
stochastic process, are reflected in the quantum channel? 
What does it imply about the true quantum nature of a 
quantum channel? What types of quantum stochastic 
processes can be simulated by classical computers? 
What are the precise obstacles for establishing a one-to- 
one correspondence between these two domains, at least 
in those cases where we have a complete characterization 
of quantum channels, like the qubit channels? In this 
letter we provide results which shed light on these 
questions. The point of view adapted in this letter, 
complements many other attempts [5| for understanding 
various aspects of open quantum system dynamics. 
Specifically we show that 



a) Four dimensional probability vectors can be charac- 
terized in complete analogy with qubit density matrices, 
where here Bloch tetrahedron replaces Bloch sphere, 

b) four dimensional stochastic matrices can be charac- 
terized in almost complete analogy with qubit channels, 

c) from these two results we make a map between four 
dimensional stochastic matrices and qubit channels, a 
map which relates many of the properties, i.e. spectrum, 
divisibility etc. of one domain to the other. 

d) We also show that the map can be extended 
to relate (P dimensional stochastic matrices to d— 
dimensional quantum channels. This map can play a 
useful role in characterization of quantum channels in 
arbitrary dimension, where a complete characterization 
is missing. 

It is worth noting that we go beyond the existing re- 
lation between d dimensional stochastic matrices and 
the restriction of quantum channels on d— dimensional 
diagonal density matrices (In these cases, the image of 
a stochastic matrix is always a zero-determinant chan- 
nel.) We show that the natural relation is between d^ 
probability vectors and c?-dimensional density matrices, 
a relation which is highly nontrivial and very reach in- 
deed. We start with a detailed exposition in the d = 2 
case. 

Consider a system with four configurations. By P^, 
/i = 0, 1, 2, 3, we denote the probability of the system 
to be in each of these configurations. The state of the sys- 
tem is described by the vector |P) = ( Pq Pi P2 P3) , 
where T stands for transpose. Obviously < < 1 and 
P^ = 1. So there are three independent parameters 

in |P). The constraint X]^=o^'m = 1 is satisfied if we 
rc-parametcrizc |P) in the following form 



1 + ri + r2 + ra \ 
1 + ri - r2 - ra 
1 - ri + r2 - ra 
V 1 - ?'i - 7-2 + 7-3 J 



(1) 



where are three real numbers, the domain of which 
will be determined later. The components of |P) can be 
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written as 

P^.^\{l + e,,.r), (2) 

where the four 3-dimensional vectors (/i ~ 0, 1, 2, 3) 

eo := (1, 1, 1) ei := (1,-1,-1) 

62 (-1,1,-1) 63 := (-1,-1,1), (3) 

point to the four corners of a tetrahedron which we call 
the Bloch Tetrahedron and denote it by A, figure ([T]). 
These vectors have the following inner products, 

e^,-e^^4:6^^^-l. (4) 
Moreover they satisfy 

XI " ^' = (5) 

From these 3-dimensional vectors we can construct an 
orthornormal basis for the four dimensional space 

K)^lil,e^f, (6) 

and write \P) as 

1^) = ^(|eo) +ri|ei) +r2|e2) +r3|e3)). (7) 

In order that \P) be a valid probability vector, all its 
entries should be non-negative and less than unity. This 
yields the following conditions for the vector r, 

- 1 < e„ • r < 3, (8) 

which means that the vector r, should lie in the Bloch 
tetrahedron A. In fact A plays the same role for classi- 
cal probability vectors as Bloch sphere plays for density 
matrices. There is a one-to-one correspondence between 
classical probability vectors and points in the Bloch tetra- 
hedron. The point r = corresponds to the completely 
mixed state, where all configurations are equally proba- 
ble. The corners of the tetrahedron correspond to pure 
states, where only one configuration has non-zero proba- 
bility. 

The basis {|e^)} allows us to characterize a four di- 
mensional stochastic matrix in an elegant way, which is 
structurally very similar to the characterization of qubit 
channels, albeit with important differences which will be 
discussed later on. A stochastic matrix Q has the prop- 
erty that the entries in each column sum up to 1, that is 
{eo\Q = (cqI, hence the expansion 

Q = \eo){eo\+U\e,){eo\+A,j\e^){ej\, (9) 

where a summation over repeated indices is understood. 
Note that Greek indices run from to 3 and the Latin in- 
dices from 1 to 3. We denote it by Qt,A, where t is a three 




FIG. 1: Bloch Tetrahedron: Every vector in this tetrahedron 
corresponds to a probability vector according to The 
corners correspond to pure states and the center to the com- 
pletely mixed state. 

dimensional vector and A is a real three-dimensional 
square matrix. When acting on a probability vector, this 
matrix induced an afiine transformation on the vector r, 
that is, 

Qt,A:|P(r))^|P(Ar + t)). (10) 

Therefore Qt,A will be a stochastic matrix provided that 
the affine transformation r — > r' :~ Ar + t maps Bloch 
tetrahedron to Bloch tetrahedron (We remind the reader 
of a similar requirement of positive maps for which 
the affine transformation should map Bloch sphere 
to Bloch sphere). This puts severe constraint on the 
matrix A and the translation t. For simplicity, hereafter 
we restrict ourselves to doubly stochastic marices, 
for which t ~ Q. In order to find the conditions on 
the matrix A, we first try to put Qo.A into a normal form. 

Normal form of Q: Consider a doubly stochastic ma- 
trix Qo,A- We can make a singular value decomposition 
of A as 

A^SAdT, (11) 

where the diagonal matrix A^i = diagonal(Ai, A2, A3) 
contains singular values of A. S and T are three dimen- 
sional orthogonal matrices and their embedding in four 
dimensions, in the basis {|e^)} is of the form 

We will then have 

Qa.A ~ SQo^AdP^ (13) 

where Qn :~ Qo.Ad i^- the right hand side is called the 
normal form of the stochastic matrix Q and has the fol- 
lowing expression 

:= |eo)(eo| + A,|e,)(e,|. (14) 

Let us now ask under what conditions, the normal 
form Qn is a stochastic matrix. This is investigated in 
the next theorem. 
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Theorem 1. A matrix Qn is a doubly stochastic matrix 
if and only if the vector X := (Ai, A2, A3) belongs to A. 



Despite the non-positivity of the matrices {A^} we 
have the foUowing 



Proof. Using ([5]), we find that the exphcit form of the 
matrix elements of Qn- They are ah of the form + 
■ A) for some e^. Therefore the elements of Qn are 
positive if and only if— l<A-e^<3, V^, which means 
that A e A. □ 

The above representation, as we will see later on, is 
very apt for constructing a map between classical four 
dimensional stochastic matrices and qubit channels. To 
this end consider the following two-dimensional matrices 



a), /i = 0,l,2,3. 



(15) 



These matrices are Hermitian, have unit trace and form a 
basis for the space of 2 x 2 matrices, but are not positive. 
Moreover they satisfy 



tr{Af,A^) = 2S^,^^, 



21. 



(16) 



Definition 1. Consider now the following map, Q 
£q, which takes a matrix Q to a qubit map 



(17) 



This map has the following properties: 
i) It respects composition, that is 

^QlQ2 = £qi^Q2^ 



(18) 



ii) It maps an orthogonal matrix like S* to a unitary 
channel, i.e. 



(19) 



where U is the representation of the orthogonal matrix 
S on the basis A^,, i.e. UA^U^^ ~ Sfj,^„A„, 

iii) If J2^,=oQ^^^^^ = ^ (E,.=o <9m,i' = 1)' t^icn £q is 
trace-preserving (unital), finally 

iv) It preserves the spectrum, that is Q\v) = v\v), if 
and only if £q(X„) = vXy, where Xy := J2a ^a^a- 

From (i), (ii) and the decomposition Q — SQnT, we 
find that 



(20) 



Note that V is defined similarly to [/, i.e. VA^V^^ = 
T^^uA,y. Equation ([20| shows the standard decomposi- 
tion of any CPT map for qubits Q. Using jlS]) and (IT71) . 
the explicit form of £q^ [Sj turns out to be 



= nitiip)! + Ktr{p(Ti)(Ti) 



(21) 



Theorem 2. The map £q^ is a unital quantum channel, 
if only if Qn is a doubly stochastic matrix. 

Proof. The Choi-Jamiolkowski Q matrix of £q_^ is found 
to be 



(22) 



The eigenvalues of this matrix can easily be calculated 
and can be expressed as 



i(l-HA-e^), M = 0,l,2,3. 



(23) 



Therefore we have arrived at the nice result that the 
eigenvalues of the Choi-Jamiolkowski matrix of £q^ are 
exactly the matrix elements of the doubly stochastic ma- 
trix Qn- Hence, these eigenvalues are non- negative if and 
only if the entries of Qn are non-negative. Therefore Qn 
is a doubly stochastic matrix, if and only if £q^ is a com- 
pletely positive unital map. □ 

Moreover in view of (pil)) we find that the map £q cor- 
responding to Q = SQnT, is also CPT and unital, even 
though for some parameters, the matrix SQnT may no 
longer be stochastic, since S and T are not necessarily 
stochastic. 

As a very simple example, the classical depolarizing 
matrix Q = p \ — p)|eo) (eo| is already in normal form 
and is mapped to the quantum depolarizing channel 

The above results suggest that we may also be able to 
map the generators of classical stochastic matrices in four 
dimensions to Lindblad generators [l^, U | of Markovian 
quantum evolution for qubits. This is indeed the case, 
provided that the matrix generates symmetric stochastic 
processes, i.e. = Qv.ii- 

Let H be such a symmetric generator which, due to the 
condition (eo|iJ = 0, will have the following expression 
in the basis ([5]) 



(24) 



where Hij = Hji . In the computational or more aptly the 
configurational basis, the off-diagonal elements of such a 
matrix, representing the rates of jumps between different 
configurations are non-negative and consequently the di- 
agonal entries are negative. It also has a normal form as 
H = S'HnS'^, where 



Hn = y^fei|ej)(ei 



(25) 



4 



and hi are the singular values of H . S' is of the same 
form as in Wc now have 

Theorem 3. The matrix Hn is a generator of a clas- 
sical stochastic process, if and only £h„ is the Lindblad 
generator. 

Proof. The off-diagonal entries of Hn are all of the form 
h-Ci (for 1 = 1,2, 3), where h := (hi, /i2, /13) and the diag- 
onal elements are all equal to /i • cq. A similar calculation 
which led to (PO]) . now leads to 

^-ffn (P) = tr{pa.i)ai. (26) 

i 

We now resort to a theorem in [l2j according to which, 
for a map p — > L{p), L is a Lindblad generator if and 
only if it is Hermitian, L*{I) = 0, where L* denotes 
the dual map, and oj'^L^uj'^ > 0. In this last expres- 
sion, uj^ = I — where |w) = "^(|00) + |11)) is a 

maximally entangled state, and L is the linearized form 
of the map, when it acts on the vectorized form of p, 
Pi_j — > Lij^kiPk.i and finally (.)'" represents the invo- 
lution {L^)ij^ki = Likji- The first two conditions are 
obviously satisfied. For the third we note from that 
= ^ hi^cTiCSial) that the eigenvalues of uj^Sj^^uj^ 
are exactly equal to ■^h- Ci, which are nothing but the off 
diagonal elements of H. So £h„ is a quantum Lindblad 
generator if and only if Hn is a classical generator. □ 

Corollary 1. For any matrix of the form H ~ S'HnS'^ , 
where S" is as in the corresponding £h will also be 

a Lindblad generator. 

Proof. Let U be as before, (i.e. after equation ((HI)). 
Then it is clear that 

£h{p)^U~^£h,SUpU'^)U, (27) 

Using the relation £^ = 2T(f ) we find 

£^h^{U®U*)-^£^hSU®U*). (28) 

Now since a;-'- commutes with U ® U* , we find 

^ r - r 

u:^£h,-, > if and only if uj'^£h > 0. Therefore 

the six-parameter family of the matrices H also lead to 
quantum Lindblad generators. In such a case, if Q = e^, 
then we find that £q = e^" , that is H and £h gener- 
ate respectively the classical stochastic process and its 
quantum image. □ 

Returning to our map Q — 5- £q , we would like to elab- 
orate on the decomposition of stochastic matrices and the 
effect it has on their quantum counterpart. Consider the 
normal form Qn- For A to lie inside the tetrahedron, the 
magnitude of components should be less than one. The 
matrix Qn shrinks the Bloch Tetrahedron in the principal 
directions by a factor A^. The transformed tetrahedron 



still lies inside A. However rotations and reflections in- 
duced by 5* and T will bring small parts of the new trans- 
formed tetrahedron to outside A. The reason is that the 
matrices 5* and T are only orthogonal and not stochas- 
tic [1^. Therefore while SQnT may be non-stochastic, 
its quantum counterpart U o £q^ o V is certainly a qubit 
channel. 

Conversely there may be stochastic matrices for which 
A does not belong to A, whose quantum image is not a 
qubit channel. As expected, there shouldn't be a one-to- 
one correspondence between classical and quantum pro- 
cesses, nevertheless we have found a map from stochastic 
matrices, whose image has a very large overlap with qubit 
channels and have also identified why such an overlap is 
not complete. It seems that two types of mapping are 
possible, one can either use orthnormal and hence non- 
positive operators, (like := ^{I + e^ ■ cr) in d = 2) or 
else choose positive and non-orthonormal operators, (like 
:= + ■ a)). We think that these two kinds 

of maps complement each other, and they are induced 
from two types of mapping the Bloch tetrahedron to the 
Bloch sphere, a subject which we will explore in more 
detail elsewhere [l^ . 

To what extent these considerations can be general- 
ized? Let {A^} and {i?^} be two basis sets of Hermitian 
operators in M^2, (the space of d^ x d^ matrices) with 
the following properties: 

^A^ = /, tr{B,) = l. (29) 

We can then define a map $, which takes any Q G 
to an £q defined as 

£q{.p)-^Y.Q>^'MpA.)B^. (30) 

If the operators {A^} and {B,^} are positive, (which will 
no longer be orthonormal), then $ maps any stochastic 
(and doubly stochastic) matrix to a CPT (and unital) 
map, readily verified by forming the Choi-Jamiolkowski 
matrix of £q. It also follows the composition rule 

^Qi^Q-i ^ £qiGQ2^ (31) 

where G is the stochastic matrix given by G^.jy = 
tr(A^B,). 

In summary we have suggested a way for mapping 
classical stochastic matrices to quantum channels in 
d dimensions. In particular, in d = 2, we have shown 
that four dimensional probability vectors and stochastic 
matrices, can be described in almost complete analogy 
with qubit density matrices and qubit quantum channels. 
This similar structure allows us to make a one-to-one 
correspondence between the normal form of doubly 
stochastic matrices in four dimensions and the normal 
form of qubit quantum channels. These findings may 
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be expanded in a few directions. Besides considering 
the case of non-doubly stochastic matrices, it obviously 
begs for a generahzation to arbitrary dimensions, a 
task which is highly non-trivial. The immediate task 
is to generalize it to the case of multi-qubits, where a 
uniformly spaced basis like © exists due to the existence 
of Hadamard matrices in 2*^ dimensions. Once this is 
done, one can also see how classical correlations lead to 
correlated quantum channels. Moreover these findings 
and the homomorphic property (|18|) . may shed light 
on the problem of divisibility and Markovianity [Hj of 
classical and quantum maps and relate this problem in 
these two domains. 
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